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Using linear kinetic plasma theory the relation between electron density and mag- 
netic field fluctuations for low-frequency plasma waves for Maxwellian background 
distribution functions of arbitrary temperatures in an uniform magnetic field is de- 
rived. By taking the non-relativistic temperature limit this ratio is calculated for 
the diffuse intercloud medium in our Galaxy. The diffuse intercloud medium is the 
dominant phase of the interstellar medium with respect to radio wave propagation, 
dispersion and rotation measure studies. The differences between the relation of 
C > ' electron density and magnetic field fluctuations from the linear kinetic theory as 

compared to the classical MHD theory are established and discussed. 

o 

1. Introduction 

o ■ 

Important input quantities for the quasilinear test-particle description of cos- 
mic ray transport in weakly-turbulent astrophysical plasmas are the wavenumber 
power spectra of magnetic field fluctuations. Within the plasma wave viewpoint, the 
plasma irregularities are usually modelled as a superposition of linear waves well be- 
low the ion cyclotron frequency, such as AlfVen and magnetosonic waves. However, 
the observed turbulence properties in the more distant interstellar and intcrgalactic 
plasmas are obtained from radio propagation measurements as dispersion measures, 
rotation measures and interstellar scintillation. These turbulence diagnostics are 
biased towards the high-density ionized interstellar phases with large volume filling 
factors, i.e. the diffuse intercloud gas and HII envelopes. In particular, dispersion 
measure and scintillation data are primarily diagnostics of density and only secon- 
darily of magnetic field. These diagnostics demonstrate the existence of interstellar 
density irregularities with Kolmogorov-type (oc lu~ s , s — const.) frequency power 
spectra extending over 11 decades in frequency, much below (lo << tt p ) the proton 
gyrofrequency (Rickett 1990| , Armstrong et al. 1995). 



Often the electromagnetic fluctuations are described within magnetohydody- 



namic (MHD) theory (e.g. Sturrock |1994[ Goldreich & Sridhar |1995| , Hollwcg 
1999| , Lithwick & Goldreich 2002) which is appropriate at large turbulence scale- 



length I > Imhd- However, the plasma parameter of the diffuse intercloud medium 
g = v ee /uip te ~ 10 -10 is much smaller than unity, so that a kinetic description of the 

t Permanent address: Department of Geological Sciences, University of South Carolina, 
Columbia, SC 29208, USA; black@geol.sc.edu 



2 



R. Schlickeiser and I. Lerche 



electromagnetic turbulence seems to be necessary. It is the purpose of the present 
paper to provide the relation between electron density and magnetic field fluctu- 
ations on the basis of the linear kinetic plasma theory. In paper II of this series 
we will use the kinetic turbulence relations to calculate frequency power spectra 
of electron density fluctuations from anisotropic power spectra of magnetic field 
fluctuations in the form of Alfven and magnetosonic waves. Such anisotropic in- 
terstellar magnetic field power spectra are required in order to be in accord with 
the h eating/cooling balance of the diffuse intercloud medium (Lerche & Schlickeiser 
200l|). 



2. Diffuse intercloud plasma parameters 

According to Rickett ( [1990D , Spangler ( |1991[ ), Lerche & Schlickeiser ( |200l| ) the 
diffuse intercloud gas can be well modelled as an isotropic Maxwcllian plasma with 
the following plasma and turbulence parameters: Bq = 1.3/iG, n e — 0.08 cm~ 3 , 
s = 5/3, v e ~ 2 • 10 7 cm s" 1 , T e = T p , v p = u e (m e /m p ) 1 / 2 = 4.7 • 10 5 cm s" 1 , 

kmin — 2"7f j ^Imax — 2 ' 10 Cm with Imax = 3 ' 10 Cm, k max — 27T / Imin — 

8 • 10~ 8 cm -1 with l m i n — 8 • 10 7 cm. l m in and l ma x refer to the smallest and 
largest scale turbulene scale length. We then infer for the Alfven speed Va = 10 6 
cm s _1 , so that the ratio \ = Va/v b = 0.05. The non-relativistic electron and 
proton gyrofrequencies are il e — 18 Hz and fl p = 0.01 Hz, respectively. For the 
thermal electron and proton gyroradii we obtain R e = i; e /|S! e | = 1.1 • 10 6 cm and 
R p = {mp/m^^Re = 43i? e = 4.7 • 10 7 cm. 

Because one is dealing with waves representing departures from a homogenous 
system, there are fluctuations in particle number density and magnetic field fluc- 
tuations that are correlated. The basic procedure to identify the se va riations has 
been originally given for the solar wind plasma by Wu & Huba ( 1975 ). However, 
their analysis deals with a specific range of parameter values, in particular they 
consider the case v e « Va opposite to ours. Here we calculate anew the relation 
between electron number density fluctuations, 5n e , and magnetic field fluctuations, 
SBi, for a more general range of plasma parameter values. 



3. Kinetic derivation of fluctuations in the electron number 
density 

3.1 General reduction for isotropic background distribution functions 

Expanding the Fourier-Laplace-transformcd Vlasov equation and Maxwell's equa- 
tions for collisionless plasmas to first order in perturbed quantities one obtains the 
fluctuations, Sf a , to the particle distribution function of plasma species a and the 
magnetic field fluctuations SB in terms of the fluctuations E in the electric field 
for any wave type any isotropic equilibrium distribution function fa°\p) (see e.g. 



Schlickeiser 2002, Eqs. (B.16) and (8.3.18)) 



Xf ,f , q a sgn(q a ) 7 dfi (p) f , . . , 

°fa{k,uj,p) = / da(p± [E x cos{4>-a)+EySm(<p-a)\ 

M a p dp J_ oosgn(ga ) 

+ E z p\\) exp[— [{k\\v\\ — uj)a — k±v±(sm((f) — a) — sin </>)]] (3.1) 

J 'a 
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and 



SB(k,uj) = -kx E 

UJ 



(3.2) 



where we have adopted a cartesian coordinate system with the wave vector k — 
(fcj.,0, fcii) = fc(sin#, 0, cos 9) and the uniform background magnetic field Bq = 
(0, 0, Bo) lying parallel to the z-axis. 

By integrating Sf a over all particle momenta p — jm a V, one obtains the fluctua- 
tions in number density, Sn a , for each particle species, thereby directly relating the 
number density fluctuations to the electric field fluctuations 



8n a (k.u) = J d 3 pSf a (k 1 uj 1 p) 
q a sgn(q a )(m a c) 2 



i-2-k 




/ dc t 


b f dBsmO [ 
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Jo Jo 



poo _ 

/ dx vT 

Jo 



+ XX 



2 2dfi°\x) 



dx 



L 



oosgn(q a ) 



isgn(g a )wVl + x 2 a , i n 

aa exp[ — J (Ji^ cos a — E y sm a J sin H cos < 

* 'a 

-\E X sin a + iJ y cos a] sin O sin + E z cosQ) 



exp [ — — — [k\\ a cos + k± ( 1 — cos a) sin 8 sin cf) + k± sin a sin 9 cos (/>)]] (3.3) 

* 'a 

where we have introduced the normalised momentum x = p/(m a c) so that 7 = 



Applying the general integral (Gradshteyn & Ryzhik 1965, Eq. (4.624)) 

r2-rr 

I cos 6 + D sin 8 sin 6) — 



/>7T /■ Z7T 

J(A,B,D)= dQ d(t>F(AcosQ + Bsinei 
Jo Jo 

2vr y dtF(Rt) = 4nsmh(R)/R 



(3.4) 



to the function F(t) — exp(t), where i? = Vj4 2 + B 2 + D 2 , and derivatives of Eq. 
(|3.4|) with respect to A, £?, and Z?, respectively, we reduce Eq. (3.3) to 



5n a = 



Amcq a (m a c) 7 
«1 



/>oo 

/ dx yl + x 2 x 
Jo 



dx 



isgn(q a )ujv / l + x 2 a u 
da exp[ — J (fe|| E z a+ 

oosgn(q a ) "a 



k±[E x sin a + ^(cosa — 1)]) 



sin Gx cos Gx 



Q3 X 3 Q2 X 2 



where 



c 



G = ^-Jk 2 a 2 + 2k\ (1 - cos a) 



(3.5) 
(3.6) 



3.2 Maxwellian background distribution functions 
Adopting the thermal equilibrium distribution 



,1 a 



(0) _ 



fJ. a n a 



4ir(m a c) 3 K 2 (p a ) 



cxp[-^ a Vl + x 2 



(3.7) 
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with the temperature parameter 

m„c 2 



(3.8) 



Eq. (3J5) becomes 



iq a f. a n a 



VL 2 a m a K 2 {p,a) J-oc Sg n(q a ) 



daG 3 (k\\E z a + kj_[E x sina + E y (cosa — 1)]) 



dx exp[— Q\J 1 + x 2 ] 



where 



Q = A«a + 



x sin Gx — Gx 2 cos Gx 
isgn(q a )uja 



(3.9) 



(3.10) 



The x-integrations can be readily performed using (Gradshteyn & Ryzhik 1965| , Eq. 
(3.914)) 



J(G) 



dx e 



^co S Gx = Q^^l 



x/Q 2 + G 2 

and its first and second derivative with respect to G yielding 



(3.H) 



<ix exp[— Q\/ 1 + x 2 ] 



; sin Gx — Gx 2 cos Gx 



9J . ^a 2 j 

Gr - 



dG dG 2 



QG 



so that 



Sn„ = 



iq a ^ a n a 



fl 2 a m a K 2 (Ha) J-oosgn(q a ) 



■ K 3 (^Q 2 + G 2 ) 
(Q 2 + G 2 f/ 2 

da(knE z a+ 



(3.12) 



k±[E x sina + ^(cosa — 1)])Q 



K^Q 2 + G 2 ) 



(3.13) 



(Q 2 + G 2 ) 3 / 2 
which is exact for all temperature values /i a . 

3.3 Nonrelativistic temperature limit \x a >> 1 

For non-relativistic temperatures fi a » 1 and low- frequencies lo << Q a << 
^ a Sl a we use the approximations K v (x) ~ y^r/2xexp(— x) and 

G 2 « til 



to derive for Eq. ( 3.13 ) the compact expression 

8n a ^ iq a 
n a n 2 rr. 



{k\\E z H 1 + k ± E x H 2 + k ± E y H 3 



(3.14) 
(3.15) 



It remains to solve the three a-integrals 



H, = 



H 2 = 



oosgn(q a ) 



daa exp 



isgn(q a )uja 



c 2 kW 



oosgn(o a ) 



da sin a exp 



n a 2n a fl 2 a 
isgn((j a )wa c 2 fc 2 a 2 



+ p a (cosa — 1) , (3.16) 



a, 



2^a^ 2 



+ Pa(cOSQ-l) (3.17) 
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and 



H 3 = 



°osgn(<j a ) 



da (cos a — 1) exp 



isgn((j a )wa c 2 k 2 a 2 



_d_ r° 

dpa J-c 



da exp 

-oosgn(<j a ) 

where we have introduced the kineticity 



2 i ' 

isgn(q a )wa c fcj 



+ p a (cos a - 1) 



2^ 2 



+ p„(cosa-l) (3.18) 



a - C2kl W 



(3.19) 



where R a = v t h, a /\^a\ denotes the gyroradius of thermal particles of species a. 
Using the identity in terms of the modified Bessel function I n (p) of the first kind, 



ePa cos« = J„(p )e» 



(3.20) 



we obtain 



m = - 



c 2 k 2 



E MPa) 



1+ It 



2c 2 fc 2 



c M^(j^l - -) 2 ] [«**(«.) + erf (m/^S(^T - n))] 



Likewise 



#2 



2z 



E J »(p.) y 



2c^kjj |^ a | 



(3.21) 



oosgn(q a ) 



da[e M - e~ M ] exp 



^ 2 fc 2 a 2 



2p Q fi a y \n a 



sgn(<? a )(cxp 



(T7T7-»-ir 



— exp 



<^-" +1 > 2 ]) 



Finally 



+ exp 



exp 



ic 2 k 2 |r^ a 



(^--n-1) 2 erf(* 



'^f (— -n-1)) 
2c^kn |^a| 



" + 1)2 l crf « ( ^- n+1)) 



(3.22) 



a(e-"-/„(p„)) r 



9p a 



oosgn(<j a ) 



da exp 



c 2 fc 2 a 2 



2p a ^ 2 N^a 



irp a nl ^ 9(e Pa I n (Pa)) r Ma^L ^ 



2c 2 fc 2 



E 



9pa 



exp 



i 2c 2 kn \£la 



(TOT"")' 
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sg » b „ )+rf(l /^(" _ n)) 



(3.23) 



3.4 Electron density fluctuations 

According to Sect. 2 we obtain for the electron kineticity (3.19) in the diffuse 
intercloud medium 



sin 2 9 , . „ 1 , . o „ 

{R e k) 2 < -{R e k max f = 2.6 • 10~ 3 « 1 



(3.24) 



3 y ^ ' - r 

values much smaller than unity at all wavenumbers and propagation angles. For 
electr on fluctuations it is therefore justified to take the limit p e — > of Eqs. ( 3.21 )- 



It is convenient to introduce Dawson's integral (Lebedev 1972 , p.l9ff. 



D\x] 



dte t 



o 



and the parameters 



and 



2 _ /i e W 



|n e |^, 



2 

3w 2 



-e x erf(ix) 



2c 2 fc 2 2v 2 k 2 



1 



2 R e kcost 



» 1 



In the limit p e — > we obtain 

ffi(p e = 0) = a 2 [2^ e C[V> e ] + wr 1 / 8 ^ - * 2 



(3.25) 

(3.26) 

(3.27) 
(3.28) 



# 2 (p e = o) = iH/, e 



DiMi-—}-D[Mi + — 



[exp(^ e 2 (l + i^) 2 ) - exp(-^ 2 (l 
2 u) 



e| \2 



) 2 )] 



and 



Z?[a e - f/'e] + £>[«e + ?/> e ] - itt 1/2 sinh.(2a e ^ e ) exp(--0 2 - a 2 

a. 



(3.29) 



^ 3 (Pe = 0) 



ir 1 / 2 in. 



' r , lOeL r , |O e 



2£# e 



2e -^ - exp(-^ 2 (l + ^) 2 ) - exp(-^ 2 (l - ^) 2 ) 



D[a e + ip e ] - D[a e - ip e ] - 2D[^ e ] 



+ 



2 



1 — cosh(2a e V'e) exp(— a 2 ) 



(3.30) 



These equations can be approximated further using the known properties of Daw- 
son's integral. Dawson's integral satisfies the linear differential equation 

dD(x) 



dx 



= 1 - 2xD(x), 



(3.31) 



has a maximum D m = 0.541 for x — xm = 0.924 and an inflection point at 
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= 1.502 where D — D w = 0.428. At large arguments (Schlickeiser & Mause 



1995) 



whereas at small arguments 

2 4 

D(x«l)~x x 3 -\ x 5 (3.33) 

v ' 3 15 y ' 

To proceed we have to calculate the parameter ip e from the dispersion relations of 

Alfven and magnetosonic waves in the diffuse intercloud medium. 

4. Dispersion relationships in the diffuse intercloud medium 

According to Sect. 2 the plasma parameters of the diffuse intercloud medium are 
in a range where 

v 2 p « Vl « v 2 e (4.1) 
We noted already (see Eq. ( [3.24 ) that the electron kineticity is much smaller than 



unity. For the proton kineticity 

p p = I% e =4.77(-^) 2 sin 2 (4.2) 

we also find values much smaller than unity if we limit our discussion to wavcnum- 
bers much smaller than k « 0A6k max . In this limit according to Sitenko ( |l967| . 
p. 115ff.) two low-frequency transverse plasma modes exist: the Alfven wave and 
the magnetosonic wave. 

The Alfven mode calculated from the dispersion relation detAy = obeys 



= VlM = Vlk 2 cos 2 6 Alfven mode (4.3) 

The Maxwell operator Ay in the relation AijEj = also specifies the wave's polar- 
isation characteristics 



E y _ AnA 23 + A 12 Ai 3 E z _ K{ 2 + AnA 



22 



(4.4) 



E x A13A22 — A12A23 E x A12A23 — A13A22 
For obliquely propagating Alfven waves Sitenko gives 

2 2 

Ea ~ E *( 1; " l r2 p tan 2 ; ~V^^ sin0cos J ( 45 ) 

at frequencies much smaller than the non-relativistic proton gyrofrequency (lu << 
Sip). 

The magnetosonic mode obeys the dispersion relation 

u 2 = Vlk 2 + [Vl + v 2 )]k 2 ± = Vl[k 2 + (1 + (i)k\] magnetosonic mode (4.6) 
if the velocity of sound is defined by means of the equalities 

2 = 2 2 Jl+2% forcos#> X J W 2 forcos0> X , , 

v * 2,P\l + ^ forcos#< X \H forcos0< X 



Note that the critical angle 6 C = arccos(x) = 87 degrees. The plasma beta is defined 
as 
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The polarisation vector for obliquely propagating magnetosonic waves is 

2 

£ M ~Ey(-% " 2 ;l;-t^-^- sing cos e) (4.9) 
\ Hp sin 9 V A Up / 

at frequencies much smaller than the non-relativistic proton gyrofrequency (lj << 
Sip). 



5. Alfvenic electron density fluctuations 

Using the Alfven wave dispersion relation ( |4.3| ) we obtain for Eq. ( 3.26| ) 



c.A 



'l.5x = 0.061 << 1 



(5.1) 



which is much smaller than unity and much smaller than the parameter a e (see Eq. 
(pi. 



In the limit ip e> A << 1 << a e Eqs. ( |3.28] )-( |3.30| ) reduce to 
Hf(p e = 0)~-al, Hf(p e = 0)~-l, tf 3 A (p e = 0)~ 



rV2 



2 <A 



so that according to Eq. (3.15) 

5n* 



II: 



{h\a\E z 



k i £L 



, tt 1 / 2 a e 
k± - . ^ 

2 Ve,A 



(5.2) 



(5.3) 



The dominant term in the bracket of Eq. (5.3) comes from the parallel electric held 
(E z ) term indicating that the electron density perturbation in oblique Alfven waves 
arises in response to their parallel electric fields. 

Making use of the polarisation properties of oblique Alfven waves (|4.5|) then yields 



8n£ 
n e 



ieE x 

Q 2 m e 



9 D 

k± — fc||Q 



rV2 



e Ujfi2 s in6>cos6> 



Q p tan 



(5.4) 



Using Faraday's induction law (3.2) we can express the magnetic field fluctuations 
(B x , B y ,B z ) in terms of the electric field fluctuations (E x , E y ,E z ) 



(B x ,B y ,B z ) = -kx (E x ,E y ,E z ) 



(5.5) 



yielding for Alfven waves 



B y = ± — (l + k 2 R 2 p )E x ~±—E x , B x 

VA VA 



■ C , m P\l/2 , r» COS0 

Z 7T xkR P —^E x , 
Va m e tan z 9 



B^~^(^xkR p ^-E x 
Va m e tan 9 



so that \B XZ \ « \B y \. Consequently, 



iV A kB y , . 
±— (sm( 



3 m v 1 



2 m e sin 9 



± 



ITT 1 / 2 VTLp COS 2 

2 m e sin ( 



(5.6) 



3 V A kB y . 

=R — (1 T' 

2ft pJ B sin(9 v 



rl/2 



(5.7) 
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The corresponding wavenumber power spectra are then related by 

y (1+ zL cos 4 g) ^ M 

n 2 40 2 sinV 9 Bp 



6. Magnetosonic electron density fluctuations 

Using the magnetosonic wave dispersion relation (4.6) we obtain for Eq. ( [3.26 ) 

ip e .M = y|x[l + (l+/?)tan 2 8} 1/2 = 0.061 \/l + 1.22 tan 2 6> 



0.061 for 9 < 42.2 

0.068 tan 9 for 5 > 42.2 

which is smaller than unity for propagation angles less than 9 C = arctan(14.78) 
86.1 degrees. 

In the limit tp eiM « 1 << a e Eqs. (|3.28[) - (|3.30[) reduce to 

Ff( Pe = 0) ~ -a 2 e , HM( Pe = 0) ~ -1, H^(p e =0) 

whereas in the limit 1 << tp e m << a e 



(6.1) 



tt 1 / 2 a e . 



2 i>, 



e,M 



(6.2) 



H^(p e = 0) 



, H™( Pe = 0) ~ -1, flf(p e = 0) 



2<M 

6.1 Small and intermediate propagation angles 9 < 



At propagation angles 9 < 9 C Eq. ( 3.15 ) then becomes 



_L" 



n e \l%m e * " 2 ^ e ,M 

Using the polarisation properties of oblique magnetosonic waves ( fi~9] 
zeE 1 ,, , , o i>„ u) . . . , to 



6n A J 



— ^-^Uku a 2 . sin 9 cos 9 + ik± 

2 m e v 11 U| O p 



O p sin 2 9 



7^2 ^ 
— 7T~ k ±- 



M 



Faraday's induction law yields for magnetosonic waves 



B , 



'k\\Ey, By 



ick cos 9 . v 



[^sin 2 ^ 



1 



sin z (9 



(6.3) 



(6.4) 



(6.5) 



B z = —k±E v 



so that \B y \ « \B X ^ Z \. Consequently. 

Sn" ieB z 



2 m e c 



v 2 , ,2 



Op sin 



rV2 



7T ' cua e 



2 



,1/ 



(6.6) 



ieB z , 3«2 cos 2 



2 m e c v 2i? 2 O p 



[1 + (l+/3) tan 2 <9] + « 



U2fc 2 



Q p tan 2 

-Bz 3 
Bo 



[l + (l+[3) tan 2 0]± 



rV2 



- <■« »s L> » - (1 + 0) sin 2 6] £ = ~ [1 + sin 2 9} |^ 



Ifiel) 



(6.7) 
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The corresponding wavenumber power spectra are then related by 
P%(k) P zz (k) 9 



n F = ^F ia+^e? (6-8) 



As an aside we note that Eqs. ( |6.6[ ) imply 

P zz (k)= tan 2 9P xx (k) (6.9) 

for magnetosonic waves. 

6.2 Large propagation angles 9 > 9 C 

Repeating the analysis at large propagation angles using Eqs. ( |6.3[ ) yields 

to? ^_B. ( 2 cos* (1 + /3Ax2 2fc2) (610) 
n e -B 3 x m e 



7. Comparison with classical MHD theory 



Eqs. (5.7), (|6.7|) and (6.10) and consequently Eqs. (5.8) and (|6.S|) are important 



modifications to the classical MHD results (e.g. Sturrock 1994, Ch. 14.1). 

Eq. (5.7) contradicts the classical MHD result Sn e — for AlfVen waves. Only 
at wavenumbers k << (2£I p /3Va) the kinetic result yields vanishing density flectu- 
ations in agreement with the classical MHD theory. 

For fast magnetosonic waves ((3 << 1) MHD theory yields 



where N = 
additional factor 



6n e ^SN _ SB Z 
+- m e n e so that n e ~ 



(3 cos 2 9 



SB Z 
B 



1 + /3sin 2 6» 
N/m,p. Our kinetic result (^ 

3, 



(7.1) 



contains the 



= -[l + 13 sin' 9] 



(7.2) 



which increases monotonically from /fci„(0) = 1.5 to /fcj n (90°) = 1.5(1 + (3). 



8. Discussion and conclusions 

Using linear kinetic plasma theory we have calculated the relation between elec- 
tron density and magnetic field fluctuations for low-frequency plasma waves for 
Maxwellian background distribution functions of arbitrary temperatures in an uni- 
form magnetic field. By taking the non-relativistic temperature limit we determined 
this ratio for the diffuse intercloud medium in our Galaxy. The diffuse intercloud 
medium is the dominant phase of the interstellar medium with respect to radio 
wave propagation, dispersion and rotation measure studies. We have found differ- 
ences between the relation of electron density and magnetic field fluctuations from 
the linear kinetic theory as compared to the classical MHD theory. 

Whereas shear Alfven waves are incompressive in MHD theory, linear kinetic 



theory yields the non-zero relation (5.7) even in the limit of vanishing electron 



kineticity p e = 0. Only at very small wavenumbers k « (2/3)|fi p | sui0/Va = 
(2/3)(ujp.i/c)sm9 the kinetic result agrees with the MHD result. 

For magnetosonic waves the kinetic ratio of the normalised density and magnetic 
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field fluctuations is modified from the MHD ratio by the factor fkin = (3/2) (1 + 
[3 sin 2 9) which is independent of wavenumber and varies within values of 1.5 and 
1.5(1 + 0). 

In the next paper of this series we will use these kinetic turbulence relations to 
calculate frequency power spectra of electron density fluctuations from anisotropic 
power spectra of magnetic field fluctuations in the form of Alfven and magnetosonic 
waves. Such anisotropic interstellar magnetic field power spectra are required in 
order to be in accord with the heating/cooling balance of the diffuse intercloud 
medium. 
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